Bright solitons in spin-orbit-coupled Bose-Einstein condensates 
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We study bright solitons in a Bose-Einstein condensate with a spin-orbit coupling that has been 
realized experimentally. Both stationary bright solitons and moving bright solitons are found. The 
stationary bright solitons are the ground states and possess well-defined spin-parity, a symmetry 
involving both spatial and spin degrees of freedom; these solitons are real valued but not positive 
definite, and the number of their nodes depends on the strength of spin-orbit coupling. For the 
moving bright solitons, their shapes are found to change with velocity due to the lack of Galilean 
invariance in the system. 

PACS numbers: 03.75.Lm, 03.75.Kk, 03.75.Mn, 71.70.Ej 



INTRODUCTION 

Solitons are one of the most interesting topics in non- 
linear systems. The most fascinating and well-known fea- 
ture of this localized wave packet is that it can propagate 
without changing its shape as a result of the balance 
between nonlinearity and dispersion The achieve- 
ment of Bose-Einstein condensation in a dilute atomic 
gas has offered a clean and parameter-controllable plat- 
form to study the properties of solitons ■ In a Bose- 
Einstein condensate (BEC), the nonlinearity originates 
from the atomic interactions and is manifested by the 
nonlinear term in the Gross-Pitaevskii equation (GPE), 
which is the mean-field description of BEC 13| . With at- 
tractive and repulsive interatomic interactions, the GPE 
can have bright and dark solitons solutions, respectively. 
Such dark and bright solitons in BECs have been studied 
extensively both theoretically (3-flojj and experimentally 

EMI. 

The developments with two-component BECs have 
further enriched the investigation of solitons in matter 
waves. The two-component BECs not only introduce 
more tunable parameters, for example, the interaction 
between the two species, but also bring in novel nonlin- 
ear structures which have no counterparts in the scalar 
BEC, such as dark-bright solitons (one component is a 
dark soliton while the other is brigh t) II Ol— 1221] , dark-dark 
solito ns [231. bright-bright solitons [2j426|, and domain 
walls j27H29j. 

Recently, in a landmark experiment, the Spielman 
group at NIST have engineered a synthetic spin-orbit 
coupling (SOC) for a BEC [3(| . In the experiment, two 
Raman laser beams are used to couple a two-component 
BEC. The momentum transfer between lasers and atoms 
leads to synthetic spin-orbit coupling 31-40|. This kind 



of spin-orbit coupling has subsequently been realized for 
neutral atoms in other laboratories 4lM44| . These exper- 
imental breakthroughs [I(J |45| have stimulated ex- 
tensive theoretical investigation of the properties of spin- 



orbit-coupled BECs 4(3- 6?|], which includes some early 
studies on solitons. For example, bright-soliton solutions 
were found analytically for spin-orbit-coupled BECs by 
neglecting the kinetic energy [68J . Dark solitons for such 
a system were studied in a one-dimensional ring [69j |. In 
this work we conduct a systematic study of bright solitons 
for a BEC with attractive interactions and the experi- 
mentally realized SOC j^HlEl. By solving the GPE 
both analytically and numerically, we find that these soli- 
tons possess a number of novel properties due to the SOC. 



In particular, we find that the stationary bright soli- 
tons that are the ground state of the system have nodes 
in their wave function. For a conventional BEC without 
SOC, its ground state must be nodeless thanks to the 
"no-node" theorem for the ground state of a bosonic sys- 
tem [7(|. Furthermore, these solitons are found to have 
well-defined spin-parity, a symmetry that involves both 
spatial and spin degrees of freedom, and can exist in sys- 
tems with SOC. 



We have also found solutions for moving bright soli- 
tons. They have the very interesting feature that their 
shapes change dramatically with increasing velocity. For 
a conventional BEC, the shape of a soliton does not 
change with velocity due to the Galilean invariance of the 
system. In other nonlinear systems, such as Korteweg-de 
Vries KdV systems, the shape of a soliton changes only 
in height and width with velocity (7lj . In stark contrast, 
bright solitons in a BEC with SOC can change shape 
dramatically from nodeless to having many nodes with 
varying velocity. The new feature arises because of the 
lack of Galilean invariance due to SOC [62j]. It is worth- 
while to note that a similar model was proposed a long 
time ago in the context of nonlinear birefringent fibers 
[72| . This shows that our results will find applications in 
nonlinear optics. 
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MODEL EQUATION 



A BEC with the experimentally realized SOC is de- 
scribed by the following GPE 



if i 



at 
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*■ (1) 



and 



where the spinor wavefunction = (^1,^2 

up-spin and ^2 

for down-spin. The nonlinear coefficient — g < is for 
attractive interatomic interactions, and we have taken 
.911 = 922 = 912 for simplicity. The SOC is realized ex- 
perimentally by two counter-propagating Raman lasers 
that couple two hyperfine ground states "Ji and ^2- The 
strength of SOC n depends on the relative incident an- 
gle of the Raman beams and can be changed 65(. The 
Rabi frequency A can be tuned easily by modifying the 
intensity of the Raman beams, a are Pauli matrices. A 
bias homogeneous magnetic field is applied along the y 
direction. We consider the case that the radial trapping 
frequency is large and, therefore, the system is effectively 
one dimensional 

For numerical simulation, we rewrite Eq. ([I]) in a di- 
mensionless form by scaling energy with HA and length 
with yJfi/mA. The dimensionless GPE is 



id t $ =[-\d 2 x 



iaa y d x 



The dimensionless parameters a 



7$ t •$]<!>. 

-K^/hA/r 



(2) 



and 



7 = N g^Jm/ (TiA) /% with N being the total number 
of atoms. The dimensionless wavefunctions $ satisfy 



fdxQ* 



2 1 
1 1 + 



$ 2 | 2 ) = 1. 



The SOC term iaa y d x 



in Eq. 



([2]) indicates that spin a y only couples the momentum in 
the x direction. The energy functional of our system is 



E= dx 



^<J> 2 | 2 



1*2 



(3) 
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STATIONARY BRIGHT SOLITONS 

We focus on the simplest stationary bright solitons, 
which are the ground states of the system. To find these 
solitons, we solve Eq. ([2]) by using an imaginary time- 
evolution method. Two typical bright solitons are shown 
in Fig. [T] One interesting feature is immediately noticed. 
There are "nodes" in these ground-state bright solitons. 
This is very different from the conventional BEC, where 
there are no nodes in this kind of ground state soliton as 
demanded by the no-node theorem for the ground state 
of a boson system. Our results confirm that this no-node 
theorem does not hold for systems with SOC [70(. 



0.15 

er 

-0.15 



- w 











o 



-0-4 • J -0.4 

8 -8 

x x 

FIG. 1: Stationary bright solitons at 7 = 1.0. The solid lines 
are numerical results and the circles are from the variational 
method. In (a) and (b), a = 1.0. In (c) and (d), a — 2.0. 



There can exist a unique symmetry for systems with 
SOC, spin-parity, which involves both spatial and spin 
degrees of freedom. The operator for spin-parity is de- 
fined as 



V - Pa,, 



(4) 



where P is the parity operator. It is easy to verify that 
our system is invariant under the action of spin-parity 
V . By direct observation, one can see that the bright 
solitons shown in Fig. [1] satisfy 



V 



<5> 2 (x\ 



(5) 



as the up component $1 has odd parity while the other 
component $2 is even. Therefore, these bright solitons 
have spin-parity — 1. In fact, all the ground state bright 
solitons that we have found have spin-parity —1. That 
the eigenvalue of V for these solitons is —1 and not 1 
can be understood in the following manner. When the 
strength of SOC a decreases to zero, the up component 
^1 shrinks to zero and only the down component sur- 
vives. Since the system becomes a conventional BEC 
without SOC, the no-node theorem demands that the 
surviving down component has even symmetry. As the 
SOC is turned up continuously and slowly, the symme- 
try of the second component should remain and the spin 
parity has to be —1. 

For a more detailed analysis of these bright solitons, 
we attempt to find an analytical approximation for the 
wavefunctions using the variational method. Motivated 
by the features of the stationary bright solitons shown in 
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FIG. 2: The relation between the number of soliton nodes 
2ir/J and the soliton width 5*. (a) Circles, squares, and stars 
are for 7 = 1.0, 1.5, 2.0, respectively. Here a increases from 
0.5 to 2.0 (the arrow direction), (b) Circles, squares, and stars 
are for a = 1.0, 1.5, 2.0. Here 7 increases from 1.0 to 4.0 (the 
arrow direction). The open circle and square correspond to 
the bright solitons in Figs. [Ha), andQJb) aQ d Figs. Die), and 
[2d), respectively. 



Fig. [TJ we propose the following trial wave functions for 
these solitons: 

, / Asin(27ra;/J)\ , , . 

*=Ucos(2WJ)J SeCh(s/5) - (6) 

The parameters A, B, J, and S are determined by min- 
imizing the energy functional in Eq. ([3]) with the con- 
straining normalization. The results of the trial wave 
functions are compared with the numerical results in Fig. 
[TJ where it can be seen that they are in excellent agree- 
ment. 

It is clear from Eq. (JS]) that the parameter 2ir/J can 
be regarded roughly as the number of nodes in the bright 



solitons, while S is for the overall width of the soliton. 
Both of them depend on the SOC strength a and the 
interaction strength 7. In Fig. [5] we have plotted the 
relation between 2tt/J and S, demonstrating how the 
number of nodes is related to the soliton width for dif- 
ferent values of a and 7. As shown in the Fig. [3J for 
solitons with the same number of nodes, they are wider 
for smaller interaction strength 7 [Fig. 12(a)] ; the solitons 
with the same width have more nodes for larger SOC 
strength a [Fig. (2Kb)] . 
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FIG. 3: The spin polarization (a z ) as a function of a and 7. 



It has been reported that there exists a quantum phase 
transition for the ground states in the spin-orbit-coupled 
system with repulsive interaction [65|. It is interesting 
to check whether such a phase transition exists for the 
case of attractive interaction. For this purpose, we have 
computed the spin polarization (<r z ) = J dx(&\ — $2) f° r 
these bright solitons and the results are plotted in Fig. 
[3] We see that for a given 7, it changes smoothly with 
the SOC strength a. For the cases of repulsive interac- 
tion and no interaction, the spin polarization is found to 
change sharply with a, indicating a quantum phase tran- 
sition [65[. The smooth behavior of Fig. [3] suggests there 
is no quantum phase transition. 



MOVING BRIGHT SOLITONS 

After the study of stationary bright solitons, we turn 
our attention to moving bright solitons. For a conven- 
tional BEC without SOC, it is straightforward to find a 
moving bright soliton from a stationary soliton: if the 
wave function $ s describes a stationary soliton, then 
exjp(ivx)& s (x — vt) is the wave function, up to a trivial 
phase, for a soliton moving at speed v. This is due to the 
invariance of the system under Galilean transformations. 

However, Galilean invariance is violated for a spin- 
orbit-coupled BEC [62] . To see this explicitly, we assume 
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moving solitons having the following form: 

$m{x, t) — $ v (x — vt, t) exp(iwa; — i—v 2 t) 1 (7) 

where & v is a localized function. Substitution of &m(%> t) 
into Eq. flJJ) yields 

id t <P v = [--dl + a.ay{id x -v) + a z -~f$l (8) 

Compared to Eq. @, this dynamical equation has an 
additional term ava y , indicating the violation of Galilean 
invariance. This violation means that it is no longer a 
trivial task to find a moving bright soliton for a BEC 
with SOC. 
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FIG. 4: Moving bright soliton profiles <&„(a;) from the numeri- 
cal calculation (solid line) and the variational method (circles) 
with Eq. 7 = 1.0. In (a) and (b), a = 1.0. In (c) and 

(d)a = 2.0. 



To find moving bright solitons, we numerically solve 
Eq. ([5]) using the imaginary time-evolution method. Two 
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FIG. 5: Dynamical evolution of a bright soliton under the 
influence of a small linear potential V(x) — rjx. a = 1.7, 
7 = 1.0, r\ = 0.001. For better comparison of shapes, the 
centers of the soliton densities have all been shifted to the 
zero point. The velocity of the soliton is labeled at the right 
of (b). 



typical moving bright solitons are shown in Fig. 0] where 
we see clearly that the shapes of moving bright solitons 
change with their velocities. As seen in Figs. 2{a) and 
B|b), when the velocity v is changed from 0.1 to 1, the 
density of the up component changes from having two 
peaks to having only one. At a larger SOC strength, 
such as a — 2 in Figs. SJc) and0|d), a small change in 
velocity leads to a dramatic change in the soliton profiles. 

Similar to the stationary soliton, these moving bright 
solitons can also be found with the variational method 
by minimizing the energy functional with the following 
trial wave function: 



sin ^ + pii cos 2^ 
cos + p 2 i sin ^ 



sech-, (9) 



with two new parameters p\ and p2- When p\ = pi= 0, 
we recover the stationary soliton in Eq. ©. The solu- 
tions obtained with the variational method are plotted 
in Fig. 2]and they agree well with the numerical results. 
It is clear from the trial wave function that the moving 
bright soliton has no well-defined spin-parity V . 

These moving bright solitons are adiabatically linked 
to the stationary bright solitons. To see this, we slowly 
accelerate the stationary bright soliton by adding a small 
linear potential in Eq. ([2]) integrating with a stationary 
bright soliton as the initial condition. The dynamical 
evolution of this soliton is shown in Fig. [5J where we 
see clearly how a stationary soliton is developed into a 
moving soliton with its shape changing constantly. Note 
that the centers of solitons in Fig. [5]have all been shifted 
to the zero point for better comparison between shapes. 

We note that there are other solitons, which also 
change their shapes with velocity, for example, dark soli- 



tons in a BEC and bright solitons in the KdV system [71 



This change is also caused by the lack of the Galilean 
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invariance in the system. For the dark soliton, the con- 
stant background provides a preferred reference frame 
and breaks the Galilean invariance. In the KdV system, 
the violation is caused by the non-quadratic linear dis- 
persion. However, in these systems, the change in shape 
with velocity is not as dramatic: there are only changes 
in the height and width of the solitons. With a spin- 
orbit-coupled BEC, the number of peaks in the solitons 
can change with a slight change in velocity. 

CONCLUSION 

We have systematically studied both stationary and 
moving bright solitons in a spin-orbit-coupled BEC. 
These bright solitons have features not present without 
spin-orbit coupling, for example, the existence of nodes 
and spin-parity in the stationary bright solitons, and the 
change in shape with velocity in the moving bright soli- 
tons. Although there are multiple peaks in the soliton 
profiles, the bright solitons that we have found are sin- 
gle solitons. It would be very interesting to seek out 
multiple-solitons solutions for this spin-orbit-coupled sys- 
tem. These bright solitons should be able to be observed 
in experiment. One can apply the same Raman laser 
setup to generate the synthetic spin-orbit coupling for 
an optical dipole-trapped 7 Li condensate where the in- 
teratomic interaction is attractive by nature. 
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Note added. Recently, we noticed a newly posted 
preprint on bright solitons in a BEC with SOC 73] . Some 
of their results overlap with ours. 
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